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Abstract 

The possibility of analyzing the node f2 of the GP-B satellite in or- 
der to measure also the Lense-Thirring effect on its orbit is examined. 
This feature is induced by the general rclativistic gravitomagnetic com- 
ponent of the Earth gravitational field. The GP-B mission has been 
launched in April 2004 and is aimed mainly to the measurement of 
the gravitomagnetic precession of four gyroscopes carried onboard at 
a claimed accuracy of 1%. of better. The aliasing effect of the solid 
Earth and ocean components of the solar K\ tidal perturbations would 
make the measurement of the Lense-Thirring effect on the orbit un- 
feasible. Indeed, the science period of the GP-B mission amounts to 
almost one year. During this time span the Lense-Thirring shift on the 
GP-B node would be 164 milliarcseconds (mas), while the tidal per- 
turbations on its node would have a period of the order of 10 3 years 
and amplitudes of the order of 10 5 mas. 

Keywords: GP-B, Lense-Thirring effect, orbit, longitude of the ascending 
node, orbital tidal perturbations 

1 Introduction 

The main scientific task of the Gravity Probe B (GP-B) mission 1 j, launched 
in April 2004, is the measurement of the general relativistic precession 2 of 
four gyroscopes carried onboard induced by the Earth gravitomagnetic field 
[3]. The claimed accuracy 1 is of the order of 1% or better. The duration of 
the experiment is almost one year. 

Measuring directly gravitomagnetism in a reliable and accurate way 
would be a very important test of fundamental physics. 

1 See, e.g., the leaflet Gravity Probe B Launch Companion downloadable from 
http://emstein.stanford.edu/index.html 



Up to now, the only attempts to measure it in the gravitational field of 
the Earth have been performed by Ciufolini and coworkers |1J [3] . They have 
tried to measure the secular Lense-Thirring effect [H] on the whole orbits 
of the existing laser-ranged geodetic LAGEOS and LAGEOS II satellites by 
analyzing time series of the combined orbital residuals of the nodes f2 of 
LAGEOS and LAGEOS II and the perigee to of LAGEOS II over time spans 
some years long. The adopted combination is not affected, by construction, 
by the first two even zonal harmonics J2 and J4 of the static part of the 
multipolar expansion of the Earth gravity potential whose induced classical 
precessions severely alias the genuine gravitomagnetic trends of £1 and 00. 
The claimed accuracy is of the order of 20-30% (Sj- For different, more 
conservative but, perhaps, more realistic estimates of the total error budget 
see jZllHlin]- However, the multi-satellite approach could turn out to be more 
fertile and fruitful in view of the improvements in the present and near future 
Earth gravity field solutions from the CHAMP [TO] and, especially, GRACE 
[TT1 IT2"] dedicated missions. Indeed, they could allow to greatly reduce the 
systematic error due to the even zonal harmonics, so to discard the perigee 
of LAGEOS II and use the nodes only. In 9 a J2— free combination which 
involves the nodes of LAGEOS and LAGEOS II only has been explicitly 
proposed: according to the recently released preliminary GRACE models, 
it would allow a measurement of the Lense-Thirring effect at a 15% level of 
accuracy (1 a). Recently, a combination involving the nodes of the geodetic 
LAGEOS, LAGEOS II and Ajisai satellites and of the radar altimeter Jason- 
1 satellite has been put forth [Ej: it would allow to reduce the error due 
to the geopotential by cancelling out the first three even zonal harmonics. 
The major problems could come from the fact that the non-gravitational 
perturbations on Jason-1 should be carefully dealt with. It must be pointed 
out that the multi-year multi-satellite J^-free approach allows to perform, 
in principle, as many analyses as one wants because of the extremely long 
lifetimes of the satellites to be used 2 and of the availability of data records 
many years long. Moreover, since, in this case, one is interested in the 
precise satellites' orbit reconstruction, such analyses can benefit from the 
improvements both in the observational techniques and in the modelling of 
the perturbing forces acting on the satellites which will become available 
in the future. So, the reliability and the accuracy of a measurement of the 
Lense-Thirring effect based on such an approach is not fixed once and for 
all but is increasing, at least up to certain level. The launch of another 

2 The lifetime of the LAGEOS satellites is of the order of 10 5 years. LAGEOS and 
LAGEOS II are in orbit since 1976 and 1992, respectively. 
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LAGEOS-type satellite like the proposed LARES El or, more recently, 
OPTIS would further enforce the reliability and precision of such a 
measurement, which could be pushed to the 1% level or, perhaps, even 
better. 

It seems, then, legitimate to ask if it would be possible to enlarge and 
enforce the significance of the GP-B gravitomagnetic test by measuring also 
the Lense-Thirring secular precession of its node. A previous treatment 
of this problem can be found in |I7] . Apparently, this possibility sounds 
very appealing because of the strictly polar orbital configuration retained 
during the science phase and of the notable accuracy reached by the most 
recent GRACE-based terrestrial gravity models. Indeed, the competing sec- 
ular nodal precessions induced by the even zonal harmonics of the Earth 
gravitational field are proportional to cosz [Hj, where i is the inclination of 
the satellite orbital plane to the Earth equator assumed as {x, y} reference 
plane. Although GP-B is much more sensitive to the higher degree mul- 
tipoles of the Earth gravity field than, e.g., the LAGEOS satellites 3 , the 
fact that iGP-B = 90.007 deg together with the high accuracy of the latest 
solutions of the Earth gravity field from GRACE should constrain the sys- 
tematic error due to the mismodelling in the even zonal harmonics of the 
geopotential to an acceptable level. For the combined impact of the depar- 
tures from nominal polar orbital configurations and of low altitudes see |18] . 
But, as we will show, this is not all the story. 



2 The impact of the static and time-varying part 
of the Earth gravitational field 

The secular Lense-Thirring precession of the longitude of the ascending node 
0, of a test mass freely falling in the gravitational field of a central spinning 
mass with proper angular momentum J is 

^ LT = c 2 a 3(l_ e 2)3/2' W 

where G is the Newtonian gravitational constant, c is the speed of light in 
vacuum and a and e are the semimajor axis and the eccentricity, respectively, 
of the orbit of the test particle. The gravitomagnetic shift of GP-B amounts 
to 164 milliarcseconds per year (mas yr -1 in the following). See Table \^iov 
the relevant parameters of the Earth-GP-B system. 

3 The semimajor axis a of GP-B amounts to 7027.4 km, while cilageos = 12270 km. 
The classical nodal precessions fall off as R e a~ l - e+3 ^' 2 \ 
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Table 1: Relevant parameters for the calculation of the solid Earth and 
ocean tidal perturbations for some selected tidal constituents and orbital 
parameters of GP-B. The degree 1 = 2 only tidal terms have been considered. 
The considered tidal constituents are the zonal (m = 0) lunar 18.6-year tide 
(055.565 in the Doodson notation) and the tesseral (m = 1) solar K\ tide. 
See |19j and the references therein for the quoted numerical values of the 
tidal parameters. For G, GM, R, p and k 2 the IERS values have been directly 
adopted [20]. The proper angular momentum of Earth has been calculated 
as J = I to, where the value of 20 for the Earth daily angular velocity 
u and of |2j for the adimensional moment of inertia I/MR 2 have been 
adopted. The orbital parameters of GP-B are those released in the leaflet 
Gravity Probe B Launch Companion. 



Symbol Description Value Units 



G 


Newtonian gravitational constant 


6.67259 x 10" s 


cm 3 g 1 s 2 


GM 


Earth GM 


3.986004418 xlO 20 


cm° s 


R 


Earth equatorial radius 


6378.13649 x 10 5 


cm 


J 


Earth proper angular momentum 


5.86 x 10 40 


2 —1 

g cm s 


P 


Ocean water density 


1.025 


g cm -3 




Tidal height of the K\ tide 


36.87012 


cm 


F 2 ° (055.565) 


Tidal height of the 18.6-year tide 


2.792 


cm 


k 2 (K l ) 


Love number for the K\ tide 


0.257 




k 2 (055.565) 


Love number for the 18.6-year tide 


0.315 






Lag angle for the K\ tide 


-18.36 


deg 


5 20 (055.565) 


Lag angle for the 18.6-year tide 


-56.29 


deg 


C+ 1 (K 1 ) 


Ocean tidal height for the K\ tide 


2.83 


cm 




Ocean hydrodynamics phase shift 


320.6 


deg 




Load Love number 


-0.3075 




a 


GP-B semimajor axis 


7027.4 x 10 5 


cm 


e 


GP-B eccentricity 


0.0014 




i 


GP-B inclination 


90.007 


deg 


n 


GP-B Keplerian mean motion 


1.0717 x 10~ 3 


s- 1 


Pn 


GP-B nodal period 


3 x 10 10 (1136.746) 


s (yr) 


n 


GP-B longitude of the ascending nc 


•de 163.26 


deg 


n LT 


GP-B nodal Lense-Thirring shift 


164 


mas yr -1 



2.1 The role of the even zonal harmonics 

In order to give a really pessimistic and conservative estimate of the system- 
atic error due to the mismodelling even zonal harmonics of the geopotential, 
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we will sum the absolute values of the mismodelled classical precessions (Hj 
induced by the whole range of the mismodelled Ji coefficients according to 
the variance matrix of the GGM01C GRACE-based model [TJ up to degree 4 
t = 48. It turns out that, at 1 a level 

^geopot = 0.8 mas yr" 1 , . . 

(50 g eopot/^LT = 5 X 1(T 3 . 1 ' 

If a root-sum-square calculation is performed by taking the square root of 
the sum of the squares of the various mismodelled classical precessions a 
relative error of 1 x 10 -3 is obtained. 

Another source of potential bias when a single orbital element is used 
is represented by the secular variations of the even zonal harmonics Jg. In 
HU it has been shown that an effective jf ~ j 2 + 0.371 j 4 + 0.079 j 6 + 
0.006Js — 0.003 Jio-.. can be introduced. Its magnitude is of the order of 
(—2.6 db 0.3) x 10~ n yr . It turns out that its mismodelled part would 
induce a secular drift of the GP-B node of 3.1 x 10~ 3 mas only in one year. 



2.2 The role of the solar K\ tide 

As we will now show in detail, the major problems for a measurement of the 
nodal Lense-Thirring shift with a polar orbital geometry ^H] comes from 
the classical time-varying nodal perturbations induced by the solid Earth 
and ocean tides. 

A given tidal constituent of frequency / induces the following solid Earth 
tidal perturbation on the node of a satellite |U 



An 



f 



" EE 7 ^E E — TP G ^P1 T k tL H T Sin limpqf , 



no? VI -e 2 sm i ~ ( ■ — \ a j — q — x ai Jp 



, (3) 

where g is the acceleration of gravity at the Earth equator, R is the Earth 
equatorial radius, Ae m is given by 




+ m)V (4) 

the quantities F^ mp {i) and G^ pq {e) are the so called inclination and eccentric- 
ity functions |23], kf^ and Hf are the Love numbers and the solid Earth 
4 For the LAGEOS satellites a calculation up to I — 20 is well adequate. 
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tidal heights, respectively, f p is the frequency of the orbital perturbation 
given by 

f p = (£- 2p)dj + (£-2p + q)M + m(0 - 9) + a, (5) 

with 

c- = ji& + (j2 - ji)s + hh + j 4 p + j 5 N' + j 6 Ps, (6) 
while the phase of the sinusoidal function is 

li mpq f = (£- 2p)uj + (£-2p + q)M + m(U - 9) + at - 6 imf . (7) 

The angular variable A4 is the satellite's mean anomaly, s,h,p,N' ,p s are 
the luni-solar mean longitudes, 9 is the Greenwich sidereal time and the 
jk, k = 1, ..6 are small integers which can assume negative, positive or null 
values. They are arranged in the so called Doodson number 

Jl (J2 + 5) (h + 5) . 0/4 + 5) (i 5 + 5) (h + 5) (8) 

by means of which each tidal constituent / is named. In it the integer 
ji classifies the tides in long period or zonal (ji = 0), diurnal or tesseral 
(ji = 1) and semidiurnal or sectorial (ji = 2). Finally, 5e m f is the phase lag 
angle of the response of the solid Earth with respect to the considered tidal 
constituent. 

The ocean tidal perturbation on the node of an Earth satellite induced 
by a tidal constituent of frequency / is |19| 



1 



=?-EEE A 



± 



i R 2 ■ ■ ^ ^ Imf 

na^V 1 - e / smz e=2m =o + 



x 



(9) 

COS ^W?/ Jf-m odd 



p=0 q=— oo 

where ± indicates the prograde and retrograde components, 

Af mf = AvGRp (±±±i\ Cf mp (10) 

p is the ocean water density, k\ is the Load Love number, Cf m ^ are the ocean 
tidal heights, 

-rLpgf = 2 p)" + {£-ip + q)M + m(0 -e)± (at - ef mf ), (11) 

and is the phase shift due to hydrodynamics of the oceans. 
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We will deal with the long-period components averaged over one orbital 
revolution (I — 2p + q = 0) because we have to compare their effects with 
the secular Lense-Thirring effect. In particular, we will focus our attention 
on the degree 1 = 2 terms. For I = 2, p runs from to 2, and so, in virtue 
of the condition t — 2p + q = 0, q assumes the values —2, 0, 2. Since the only 
non vanishing eccentricity function of degree I = 2 is that for p = 1, q = 0, 
it follows that the condition t — 2p = is also fulfilled. The frequencies of 
the perturbations are, in this case 

/ p = f/ + mO, (12) 

with 

f / = C?2 - m)s + j 3 h + ]4p + j 5 N' + j G p s . (13) 

Since the science phase of GP-B will be almost one year long, it is of the 
utmost importance to see if there are some tidal perturbations, with rele- 
vant amplitudes, which could resemble as superimposed linear trends on the 
gravitomagnetic signal during such a time span. The most powerful tidal 
constituents in affecting satellite orbits are the zonal (m = 0) lunar 18.6- 
year tide, with Doodson number (055.565), and the tesseral m = 1 solar K\ 
tide with Doodson number (165.555). For such tidal lines the inclination 
factor of the perturbing amplitudes are 



1 dFwx 
sin % di 

1 dF 211 
sin i di 




(14) 

sini , (15) 



respectively. This implies that a nearly polar orbital geometry will be af- 
fected mainly by the K\ tide than by the 18.6-year tide. This is an im- 
portant feature because the period of the nodal perturbation induced by 
(055.565) depends only on the luni-solar variables and amounts to 18.6 
years. The situation is quite different for the K\ tide. Indeed, from eas. ()12j) - 
(|13|) it follows that its nodal perturbation has the same period of the node 
of the satellite. For GP-B it amounts to 1136.746 years. From Table ^ 
and by assuming a time span of one year eq.(j3J) yields a nominal shift 
of AQiK^f = 2.01264 x 10 5 mas. By assuming a 0.5% uncertainty 
in the Love number Qlj one gets 5[Afi(A'i)g 1 2 d ] = 1000.6 mas. Eq.© 
yields a shift of AQ(i^i)^ff 2 an = 4.94785 x 10 5 mas after one year. By 
assuming an uncertainty of 3.8% in the ocean tidal height , one gets 
8[Att(Ki) e ^ n ] = 1.8801 x 10 4 mas. The 18.6-year tide would not pose 



7 



particular problems. Indeed, the nominal amplitude of its nodal perturba- 
tion would amount to 2.7 mas only; the error in the Love number at that 
frequency is estimated to be ^F3\ 1.5%, so that <5[Ar2(055.565)] = 0.04 mas. 
These results clearly show that the bias induced by the K\ tide on the node 
of GP-B would not allow to use it in order to measure the Lense-Thirring 
effect. Such results lead to the same conclusions of |17j . 

2.3 A multi-satellite approach 

At this point one could ask if it would be possible to include the node of 
GP-B in some multi-satellite J^-free combination. The answer, also in this 
case, is negative because of the polar geometry of its orbit. As shown in 
|18| . the coefficients with which the residuals of the orbital elements of a 
polar satellite would enter some combinations tends to diverge for i ~ 90 
deg. In the case of GP-B one could consider, e.g., a 5f2 LAGEOS + c\5(l GP ~ B 
combination. The coefficient c\ of GP-B would be equal to 5 -398 so that the 
impact of the K\ tide would be further enhanced. 

3 Conclusions 

In this paper we have investigated the possibility of measuring the gravit- 
omagnetic Lense-Thirring effect on the orbit of a test particle by analyz- 
ing also the nodal rate of the GP-B satellite in the gravitational field of 
the Earth. It has been launched in April 2004 and its main task is the 
measurement of the gravitomagnetic precession of four gyroscopes carried 
onboard at a claimed accuracy of 1%. The main problems come from the 
fact that the orbital perturbations induced by the solid Earth and ocean 
components of the solar K\ tide on the GP-B node would resemble aliasing 
linear trends superimposed on the genuine relativistic linear signal. Indeed, 
the observational time span would be of almost one year while the period 
of such perturbations is of the order of 10 3 years. Moreover, their ampli- 
tudes are three orders of magnitude larger than the Lense-Thirring effect. 
This drawback would be further enhanced by including the GP-B node in 
some multi-satellite Ji— free combinations because the coefficient with which 
it would enter them would be very large due to the polar geometry of the 
satellite. 

5 The systematic error due to the mismodelling in the even zonal harmonics would 
amount to 50-100%, according to GGM01C model. 
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